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Using an algebraic formulation, we explore two well-known degeneracies involving the 
mass-squared differences for three-neutrino oscillations assuming CP symmetry is con- 
served. For vacuum oscillation, we derive the expression for the mixing angles that permit 
invariance under the interchange of two mass-squared differences. This symmetry is most 
easily expressed in terms of an ascending mass order. This can be used to reduce the 
parameter space by one half in the absence of the MSW effect. For oscillations in matter, 
we derive within our formalism the known approximate degeneracy between the stan- 
dard and inverted mass hierarchies in the limit of vanishing S13. This is done with a 
mass ordering that permits the map A31 ^ — A31. Our techniques allow us to translate 
mixing angles in this mass order convention into their values for the ascending order 
convention. Using this dictionary, we demonstrate that the vacuum symmetry and the 
approximate symmetry invoked for oscillations in matter are distinctly different. 
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1. Introduction 

In general, neutrino oscillations are achieved by relating, via a nontrivial unitary 
mixing matrix [7, flavor states to mass eigenstates of a Hamiltonian H\ that is, 
Vf = Uvm- -^s the flavor states are not eigenstates of the Hamiltonian, flavor oscil- 
lation occurs as the particles propagate through space. The degree of mixing among 
the mass eigenstates is characterized by the mixing matrix U ; and the frequency 
of the oscillations is characterized by the differences in the squared masses of these 
eigenstates. Phenomenologists seek to determine the number of relevant neutrinos, 
the mixing matrix, and the mass-squared differences. In a two neutrino theory, one 
has one mass-squared difference, and physical arguments reduce the mixing matrix 
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to an element of the commutative group U{1). As there are only two parameters 
of interest, symmetries of the relevant formulae are readily apparent even if one 
considers the additional complications brought on by matter effects. For three neu- 
trinos, the situation is more interesting but more complicated. The mixing matrix 
is now an element of a noncommutative group S'i7(3), parameterized by three real 
mixing angles and a Dirac phase; and one now must deal with two independent 
mass-squared differences. Adding matter effects into the description, further com- 
plicates the issue. Though analytic expressions do exist in this general setting^, any 
potential symmetries or degeneracies among the parameters are mired in awkward 
relations and transcendental equations. For iV > 3 neutrinos, using traditionally 
formulated oscillation equations to determine symmetries is nearly hopeless even 
without matter effects. 

In order to determine such degeneracies, one should consider neutrino oscillations 
in the context of the algebraic formulation developed in Refs. El and El valid for a 
fixed arbitrary number of flavors. We will restrict our discussion to a three neutrino 
theory; however, the methodology is readily generalized. In Ref. |31 we utilized this 
formalism to derive symmetry relations among the three mixing angles and CP 
phase in the mixing matrix. Here we shift our attention to symmetries involving the 
mass-squared differences. 

We begin with a discussion of the discrete CP and T symmetries for oscillation in 
vacuo. The differences that we shall find between the vacuum and matter oscillations 
are rooted in these symmetries. Assuming CP invariance, we derive a degeneracy 
with respect to the mass-squared differences for three- neutrino vacuum oscillations. 
This was observed by Ahluwalial^. The symmetry is exact allowing one to reduce 
the physically relevant parameter space by one half, as in the two-neutrino case. 
In order to discuss the approximate symmetry for neutrino oscillations in matter, 
one must appeal to phenomenology for guidance. Data demonstrate that there are 
two distinct scales for the mass-squared differences with one mixing angle near zero 
1^13. We discuss this well-known approximate degeneracy within our formalism. The 
two separate degeneracies are most easily accessed through different mass ordering 
conventions. We develop a dictionary so that we may relate the two mass ordering 
conventions which allows us to demonstrate that these are two distinct symmetries. 

2. Discrete symmetries 

We refer to Ref. E] for a field-theoretic treatment of vacuum neutrino oscillations. 
The probability that an a-fiavor neutrino will be detected as a /3-fiavor neutrino at 
time t after its creation is 



The operators in this expression are defined as follows. The projection operator 
onto the flavor a is P", with components in the mass basis {P°')jk = U*jUak- Time 



(1) 
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dependence of the projection P°'{t) is as in the Heisenberg picture, 

= e'^'P'^e-'"', (2) 

where, in our case, we take to be a Dirac-hke Hamiltonian. The remaining op- 
erator P+ (or P-) projects onto the positive-energy particle (negative-energy an- 
tiparticle) subspace, 

P± = i(l±Hf-i), (3) 

where we define £ to be the positive square root of the squared Hamihonian. 

The abihty to swap the mass-squared differences is underpinned by the funda- 
mental symmetries of time reversal and CP conjugation. We begin by explicating 
how these symmetries are implemented in the formalism. As our Hamiltonian is 
Dirac, there exist two antilinear involutions of import: time reversal, denoted by 
T, and CP conjugation, denoted by S. Time reversal T commutes with the Hamil- 
tonian, and S anticommutes with the Hamiltonian. As T is an involution and the 
trace is cyclic, we have from Eq. Q 

Vo.^p{t, U) = ]^ty[TP+P'-[t)P+P^T] . (4) 
From the antilinearity of T, we achieve the equality 

rc.->p{t,u)^Vo.^p{~t,u*), (5) 

where the dependence on U is implicit in the definition of P" as noted above. 
The effect of the CP conjugation operator on the (anti-)particle projection is 

SP±I] = P^ , (6) 

due to the anticommutativity of S with H . We also have 

Ee*^*!] = e'""- , (7) 

as E is antilinear. The consequence to the oscillation probability is 

Vo.^p{t, U) = itr[EP+P"(t)P+P'3q 

= V^^p{t,U*). (9) 

From Eq. jSJ and Eq. lO, we note that the theory conserves CP and T separately 
ii U = U* . Considered jointly, we see that successive application of these operators 
indicates CPT invariance in general 

ro.Mt,u)^v^^p{-t,u). (10) 



In what follows, we will restrict our study to a CP invariant theory; i.e., we set 
U = U*. 
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3. Vacuum oscillation 

In the ultra-relativistic limit {E ^ m), the oscillation probability Eq. |^ in vacuum 
for source-to-detector distance L and neutrino energy E becomes 

Va^^piL/E) - tr[e'■^^/2Bpag-^A1L/2i^p/3J 

where Ad, in the mass basis, is the diagonal matrix with entries (m^, m|, to§). In 
this formulation, it is clear that adding multiples of the identity to A4 leaves the 
trace invariant; in other words, it is only the mass-squared differences which are 
of dynamical consequence. In this limit, time reversal invariance is equivalent to 
invariance under the map A4 i— > —M. Making the summation explicit, the trace in 
Eq. Hll|) is seen to yield the usual expression for the oscillation probability, 

3 

Pa,^p{L/E) = - 4 ^ Ua^jUa^kUpkUpj sin^i^jk) (12) 

j>k 
j,k = l 

where (fjk ■= AjkL/AE with Ajk := m| — to|. Here, the dependence on the mass- 
squared differences is readily apparent. We remark that by definition they must 
satisfy A12+A23-I-A31 = 0; hence, only two of the three are independent parameters. 

In the context of vacuum oscillation, an exact degeneracy exists between these 
two independent mass-squared differences. This degeneracy is valid for all values 
of mass-squared differences and all mixing angles; additionally, it is independent of 
the well-known approximate degeneracy for oscillations in matter. This symmetry 
is most easily accessed if we insist upon the mass ordering 

mi < 1712 < m^. (13) 

Such an ordering is entirely general as permutations of the mass labels can be 
achieved via unitary transformations which preserve the underlying physics. Equiv- 
alently, we may characterize the oscillation dynamics with the positive quantities 
A21, A32 > 0. We make no assumptions regarding the relative values of these two 
mass-squared differences, as our result holds in general. Derived herein, we may in- 
terchange the values of these two mass-squared differences and maintain equivalent 
oscillation probabilities given a suitable change in the mixing matrix U . As such, 
we need to adopt a parameterization of the mixing matrix; we use the standard 
representation ^ 

(C12C13 S12C13 S13 \ 

-S12C23 - C12S23S13 C12C23 - S12S23S13 S23C13 (14) 
S12S23 - C12C23S13 -C12S23 - S12C23S13 C23C13 / 

where Cjk = cosOjk, Sjk — sin9jk, and the ranges on the mixing angles 9jk are 
discussed in Refs. l^landUl For our purposes, it will be useful to write the mixing 
matrix as three separate rotations 

1/(^23, 013, 012) = f/l(023)f/2(ei3)f/3(012) ■ (15) 
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Here, Uj{9) represents a proper rotation by angle 9 about the jth axis in R"^. 

We are interested in the transformations of the mass-squared matrix which re- 
label the masses. These are rotations about the jth axis through an angle of n/2. 
The rotation Uj{Ti/2) leaves the jth axis unchanged and interchanges the other two 
axes. Conjugating M. with U2{tt/2), for instance, yields 

M' := U2{tt/2)MU2{ti/2)^ = diag(m2, m^, ml) . (16) 

There are five nontrivial ways to permute the three masses; however, any of these 
will disrupt the mass ordering adopted in Eq. (|13|l . In fact, we see that in the above 
example Eq. (|16() the mass ordering is reversed, m'l > m'j > mg . 

We may invoke time reversal invariance, or invariance under the map M >— > 
— A^. This is simply a restatement that oscillation probabilities depend only on 
the absolute value of the mass-squared differences \^jk\, as is manifest in Eq. H12|l . 
Given ascending masses as in Eq. ()13|l . the negative squared masses will satisfy the 
reverse inequality. As such, if we effect the flip of squared masses as in Eq. (I16II and 
then change their sign, we have a new mass-squared matrix A4' with A'21 — A32 
and A32 = A21, which is what we desire. 

The only remaining task is to determine what new mixing angles will result 
in an equivalent theory for mass-squared matrices A4 and A4' . Using the relation 
between the two in Eq. we have 

UMU^ ^UU2{TT/2y M'U2{tt/2)U'' . (17) 

Clearly, the mixing matrix in question is [/ = UU2{t^ /2)'^ . To make this of practical 
use, we must determine the new mixing angles O'^f^ in the chosen parameterization 
of Eq. (HH). 

In general, given the (real) matrix elements Ujk of some 3x3 unitary matrix, 
we may parameterize the matrix in terms of the mixing angles Ojk as in Eq. (|14|) 
by solving the following transcendental equations 

tan0i2 = [/i2/[/ii, sin6'i3 = C/13, tan 6*23 = t/23/f/33 (18) 

assuming [/ii,C/33 7^ 0. Should these matrix elements vanish, then we must deal 
with them independently. We first present these two special cases below and follow 
with the general situation. 

Case (i): sin0i3 = 0. Since data analyses yield 6'i3 small ^ ^, this is a particularly 
interesting case. From the symmetry relations for the mixing angles in Ref. H) it 
suffices to consider only ^13 = 0. The new mixing matrix may be expressed as 

u ^iJi{e23)Mei2)U2{TT/2)^ 

= f/i(^23 + 7r/2)C/2(0i2 - V2)t/3(vr/2); (19) 
hence, the map is given by 

^23-^23+^2, 0'i3 - 012 - vr/2, 0'i2 = 'r/2. (20) 
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Case (ii): cosf?i2 = 0. According to Ref.0I we need only consider 612 = 7r/2. We 
find 

U = f7i(023)f/2(f?i3)C/3(7r/2)C/2(7r/2)t 

= 1/1(^23 + 7r/2)C/2(0)t/3(V2 - 013), (21) 

so that 

0^3 = ^23+^2, ^'13 = 0, 0;2=V2-0i3. 

Case (iii): The remaining situations. We have 

U = C/l(023)t/2(013)C/3(012)^2(7r/2)t 

= C/i(023 + a)C/2(/3)C/3(7), 
where the angles satisfy the following set of transcendental equations 

tana — tan0i2 csc0i3, 
sin /3 = ~ cos 9i2 cos 6*13 , 
tan 7 = sin 612 cot 013 . (24) 

The new angles are thus 

023 =023 + a, 0'i3=/3, 0'i2=7. (25) 

One consequence of this symmetry is that, in the context of vacuum oscillations, 
we may reduce the parameter space as with the so-called "dark side" in a two- 
neutrino theory. One may choose to reduce the parameter space by placing further 
constraints on the mass-squared differences, e.g., A21 < A32. Alternatively, one may 
look at all possible mass-squared differences while further restricting the bounds on 
the mixing angles. 

4. Relating two mass order conventions 

This symmetry was readily realized for the convention of ascending masses; however, 
this choice is not the dominant one found in the literature. Though our choice was 
made without loss of generality, we demonstrate the relation to the mass ordering 
used extensively in three-neutrino phenomenology. In order to account for the over- 
whelming majority of neutrino oscillation data, the two mass-squared differences 
need to differ by several orders of magnitude. Typically, one uses the convention 

A21=Aq>0, IA31I = Aatm >0, (26) 

with Aq ^ Aatm, cf. Refs.|n|and|3 Though this convention establishes an ordering 
on the relative absolute magnitudes of the mass-squared differences, it is entirely 
general in that it allows A31 to be negative. For A31 positive, one has what is termed 
the standard hierarchy. This obeys the ascending mass order that we previously 



(22) 



(23) 
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adopted. When A31 is negative, then one has an inverted hierarchy. Such can be 
achieved with the mass ordering 

ma < mi < m2. (27) 

Although this convention has some distinct advantages, it can be conceptuaUy dif- 
ficult with which to work For the natural hierarchy, the ordering of the mass 
eigenstates is by ascending mass. For the inverted hierarchy, the ordering is given 
by Eq. (|27l) , which constitutes a change of basis for one half of the problem. 

For completeness sake, we derive the map between this formulation of the in- 
verted hierarchy and the physically equivalent one with ascending masses. Suppose 
we have mixing angles 9jk and mass-squared differences A21 = Aq and A31 = 
— Aatm- Equivalent oscillation probabilities can be got with the new mass-squared 
differences Aj^ = Aatm and A32 — A© (which obey the ordering m'/ < mj < mg 
and A32 <^ A21) and mixing angles given by 

^23 — 023 — arctan(l/ tan 6*12 sin 013), 
013= arcsin(sin 6*12 cos 6*13), 

= arctan(cos0i2/tan0i3), (28) 

whenever the operations are defined. Should sin0i2 = 0, then the new mixing angles 
are 

^23 = ^23 - 7r/2, 0i'3 = O, 01'2= 013-^2. (29) 
Finally, the case of sin6'i3 — 0, we have 

^23= ^23+^2, 0'/3 = 0i2, e;'2 = 7r/2. (30) 

5. Oscillation in matter 

The inclusion of matter effects destroys the above symmetry. We shall give a brief 
exposition as to why this is the case and then examine another well-known, but 
approximate, mass hierarchy degeneracy for the oscillation probability. 

For oscillations in matter, we need to modify the Hamiltonian to acc oun t for the 
charged-current interaction between the electron and electron-neutrino ^I^J Xo this 
end, one introduces an effective potential that operates only on the electron flavor. 
At the level of the particle interaction, CPT is taken to be conserved; however, as 
there exists a natural disparity between particle and antiparticle content of ordinary 
matter, one admits the possibility of extrinsic CPT violation as discussed extensively 
in Ref. El and the references therein. We continue to assume an intrinsically CP 
invariant theory; however, the MSW effect will lead to an extrinsic violation of this 
symmetry due to the change in sign on the effective potential for antineutrinos. 
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Oscillation probabilities for an oriented path F through the matter are given by^ 

P„^;3(F,^) = tr 



expji j Hdx^P°' 
exp|-i j Hdx^pf^ 



(31) 



where, in the ultra-rclativistic limit and modulo a multiple of the identity, the 
argument of the path ordered exponentials in this trace can be written as 

j Hdx ^ ML/2E + J A{x)dxP''. (32) 

For antineutrinos, one needs to change the sign on the potential A{x). This leads 
to an extrinsic violation of CP symmetry even though we have U = U*. As such, 
we are unable to reverse the ordering of the masses via the map M t-^ —A4, a, 
crucial step in the previous derivation. However, as is known, there exists another 
approximate symmetry between the regular and inverted mass hierarchies. 

The presence of the effective potential complicates our cause, but not entirely. If 
there exists some nontrivial subspace on which the MSW potential commutes with 
A4, then we say that the MSW effect has decoupled on that subspace. Physically, 
the MSW effect does seem to decouple, at least, approximately. Let us write the 
mass-squared matrix as A4 = diag(0, A21, A31). In the ideal limit of vanishing ^13, 
one notes that the image of the electron flavor projection lies wholly within 
the subspace spanned by the toi and m2 mass eigenstates. We denote by P3 the 
projection onto the eigenstates of mass 7713. The MSW potential commutes with 
this projection. Using orthogonal projections, we separate the integral in Eq. (|32(l 
into two commuting terms 



Hdx = (1 - P3) / Hdx + P3 / Hdx; (33) 
r Jr Jr 

to condense notation, we shall call the first operator on the RHS B and the second 
C. The second term is decoupled from the matter potential 

C^PsJ Hdx = P3ML/2E (34) 

where P3A4 — diag(0, 0, A31). As B and C commute, the following applies exp{i(i3+ 
C)} = exp{iB) exp{iC). In this limit, we consider the electron-electron oscillation 
probability in matter; that is, we set a = = e in Eq. H31(l . Clearly, exp(iC) 
commutes with the flavor observable P'^ so that the operators involving A31 cancel 
each other, making the probability independent of this parameter. 

It is this independence that prompts us to consider negative values for A31, while 
leaving unchanged A21 so as to not disrupt the MSW effect. As shown, the change 
of sign of A31 does not affect the Pf, — Vg oscillation probability in matter but it does 
yield the inverted hierarchy. We discuss the ramifications for the other probabilities, 
in particular, the vacuum oscillations. The first consideration is the effect upon the 
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remaining mass-squared difference. Under the map, we have A32 1— > — A31 — A21. In 
the relevant limit of A21 ^ IA31I, one has to a good approximation A32 1— > — A32. 
The vacuum oscillation formula, Eq. (|12() . demonstrates that the oscillations will 
remain unchanged upon these two changes in sign. As such, we have an approximate 
symmetry provided ^13 is near zero and A21 <C IA31I. 

We can use the results in Eq. (|28() to determine what this approximate symmetry 
looks like in the convention of ascending masses. Though they are of intrinsic inter- 
est, these relations allow us to compare the near degeneracy between the standard 
and inverted mass hierarchy in matter with the exact vacuum symmetry relating the 
two hierarchies. For ease of comparison, we restate the approximate symmetry for 
the inverted mass hierarchy with matter effects in the ascending mass convention, 

^23= ^23+^2, ^13 = ^12, e;'2 = 7r/2. (35) 

Likewise, the vacuum symmetry for Oi^ = as discussed in Case (i), Eq. if^ . is 

^23 = ^23+^/2, 0'i3 - 012 - 7r/2, 9[^^7r/2. (36) 

Putting aside the correction to the mass-squared difference A32 , we note the differ- 
ence between the mixing angles O"^ and d'l^ for the two degeneracies. This clearly 
demonstrates that the approximate degeneracy between the two hierarchies for neu- 
trinos propagating in matter is not a mere perturbation of the exact vacuum sym- 
metry. 

W e examine the matter mass hierarchy for a realistic case. We use a model ^1 
to fit oscillation parameters to data. The results yield a value of 613 = 0.13, a small 
but not insignificant number. The matter hierarchy has so far been derived here 
neglecting terms of order A21/A31. We define the symmetry so as to preserve the 
spacing between the levels, 

A'^, = -(A31 - A21) . (37) 

A numerical calculation yields that with this definition, the resulting oscillation 
parameters remain unchanged to the order of several tenths of a percent, even 
though 013 is 0.13. A numerical refit of the data yields physical results that are even 
more similar between the hierarchies. 

6. Conclusion 

We have examined, within our algebraic formalism, well-known symmetries which 
relate to the mass hierarchy question. We work in two different conventions for the 
mass ordering for the inverted hierarchy and provide the dictionary that relates the 
two conventions. We first investigate a symmetry which, in the absence of matter 
effects, is exact for all values of the parameters. We find this symmetry reduces the 
allowed parameter space by one-half in the absence of matter effects, just as occurs 
in the two neutrino case. We also derive the mass hierarchy symmetry which holds 
in the presence of the MSW effect in the limit of small ^13 and two different scales 
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for the mass-squared differences. Utilizing our dictionary that translates between 
the two conventions for the mass ordering, we show that these two symmetries are, 
in fact, different. For vanishing ^13, the relations for the mixing angles effecting the 
vacuum symmetry are given in Eq. (|20|l while the relations for the approximate 
symmetry in matter are in Eq. (|30|l . Finally, we find, like others that even for 
modest values of ^13 the approximate mass hierarchy symmetry is quite accurate. 
Given that the present ability to perform neutrino oscillation experiments yields er- 
rors substantially larger than the difference between the two hierarchy solutions, one 
tends to think of them as equivalent solutions when phenomcnologically extracting 
mixing angles and mass squared differences from data. The differences between the 
mass hierarchies are important in the context of looking for CP or CPT violation, 
themselves very small effects. However, there is other physics which in principle dis- 
tinguishes between thes e solutions. An example would be gravitationally induced 
neutrino oscillations or oscillation in high density systems such as supernovae. 
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